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ABSTRACT 


This  paper  deals  with  the  diffraction  of  time-harmonic 
electromagnetic  radiation  by  perfectly  conducting  obstacles 
immersed  in  an  inhomogeneous,  anisotropic,  conducting 
medium.  A  mathematical  fcarmulatlon  of  the  problem  is 
presented  which  is  applicable  to  obstacles  of  arbitrary  shape 
and  to  a  very  general  class  of  media,  and  the  existence, 
uniqueness  and  continuous  dependence  on  the  data  of  the 
solution  is  demonstrated. 


THE  STEADY-STATE  DIFFRACTION  OF  ELECTROMAGNETIC  RADIATION  BY  AN 
OBSTACLE  IN  AN  INHOMOGENEOUS  ANISOTROPIC  CONDUCTING  MEDIUM 

Calvin  H.  Wilcox 

§1.  Introduction.  This  paper  deals  with  the  diffraction  of  tlme-harmonlc 
electromagnetic  radiation  by  perfectly  conducting  obstacles  Immersed  In  an 
Inhomogeneous,  anisotropic,  conducting  medium.  A  mathematical  formulation 
of  the  problem  Is  presented  which  is  applicable  to  obstacles  of  arbitrary  shape 
and  to  a  vary  general  class  of  media,  and  the  existence,  uniqueness  and 
continuous  dependence  on  the  data  of  the  solution  Is  demonstrated. 

Electromagnetic  fields  are  represented  by  pairs  of  vector  fields,  E(x,  t) 
(the  electric  field)  and  H(x,  t)  (the  magnetic  field),  which  are  described  here 
by  their  (real-valued)  components  £j(x,  t),  Hj(x,  t)  (J  =  1,  2,  3}  relative  to 
a  fixed  rectangular  coordinate  system.  The  symbol  x  =  (x^,  x^,  x^)  denotes  a 

3 

point  in  Euclidean  space  R  and  t  Is  a  time  coordinate.  Time-harmonic 
electromagnetic  fields  have  the  form 

(1.1)  E^(x,  t)  =  Re  {e"^"‘Ej(x)},  Hj(x,  t)  =  Re  {e"^^  Hj(x)} 

where  u  is  a  real  frequency  and 

Ej(x)  «  eJ(x)  +  lEj^x),  Hj(x)  =  hJ(x)  +  IHj^x) 
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are  the  complex-valued  components  of  vector  fields  £(x},  H(x)  which  are 
Independent  of  t  .  Maxwell's  equations  for  tlme-harmonlc  fields  In  an 
Inhomogeneous,  anisotropic,  conducting  medium  filling  a  domain  (2  C  R  have 
the  form 


(1.2) 


(VXH)j-(l.Cj,+  crj^}E^  =  Ij  ' 

(VX  E)^  + 


inn  . 


Hare  the  summation  convention  is  used  (repeated  indices  are  summed  from  1  to 
3),  (VXH)j  =  aHj/ax^  -  aH^/ax^,  etc.  The  functions  ® > 

'‘jk  =  '‘jkW  and  jjj  =  O'  jjj(x)  are  real-valued  and  represent  the  components  of 
the  dielectric,  magnetic  permeability,  and  electric  conductivity  tensors, 
respectively.  and  are  complex-valued  and  represent  the  electric  and 
magnetic  current  densities. 

The  quadratic  form 

Ej(x,  t)  Ej^(x,t)  +  (ijj^(x)  Hj(x,  t)  Hj^(x,t)) 


defines  the  energy  density  for  (real-valued)  solutions  of  Maxwell  .,  equations 

in  an  inhomogeneous,  anisotropic  medium.  Hence,  the  tensors  e  „  and  u  . 

jK 

are  assumed  to  be  positive  definite.  The  conductivity  tensor  <t  ,  must  also  be 
positive  definite  if  the  medium  is  to  be  dissipative  (energy  is  absorbed,  rather 
than  created,  in  it). 
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A  perfect  conductor  (of  electricity)  Is  characterized  by  the  property  that 
the  tangential  component  of  the  electric  field  vanishes  on  Its  surface.  Thus, 

If  a  perfectly  conducting  obstacle  oC  Is  Immersed  In  a  medium  occupying 
0  =  -  O  then 

(1.3)  NX  E  =  0  on  an  , 

where  N  Is  a  normal  vector  on  an  . 

The  diffraction  problem  considered  here  (the  steady-state  diffraction  problem) 
asks  for  the  tlme-harmonlc  fields  (1. 1)  generated  by  presorlbed  tlme-harmonlo 
electric  and  magnetic  current  densities 

J(x,t)  =  Re  (e"^^  J(x)},  K(x,t)  ■  Re  {e‘^‘  K(x)> 

acting  In  the  presence  of  a  prescribed  obstacle  O  .  Maxwell's  equations  (1*  2) 
and  the  boundary  condition  (1.3)  are  necessary  conditions  on  the  solution. 

However,  they  do  not,  in  general,  determine  the  solution  uniquely.  Indeed,  if 
the  medium  Is  homogeneous  and  isotropic  and  an  has  sharp  edges  It  Is  known 
that  (1. 3)  must  be  supplemented  by  an  "edge  condition"  to  obtain  uniqueness 
[1],  Moreover,  If  n  Is  unbounded  a  "condition  at  infinity"  is  needed  to  obtain 
uniqueness  [5,  7] .  A  complete  formulation  of  the  steady-state  dt^action 
prc^lem  is  given  below  In  §2  .  The  discussion  in  the  remainder  of  this  section 
is  intended  to  motivate  the  final  formulation. 

The  time-average  of  the  energy  density  for  a  time-harmonic  electromagnetic 


field  is 


-4“ 
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Is  the  complex  conjugate  of  etc.  Most  edge  conditions 
have  been  based  on  the  physical  principle  that  the  energy  in  bounded  portion 
of  space  should  finite;  1.  e. , 

(dX  =  dXjdXjdX,) 

for  each  bounded  set  CC  R  •  This  condition  eliminates  the  possibility 

that  point-or  line-sources  of  energy  might  reside  in  a  sharp  edge.  (These 
would  clearly  lead  to  non-uniqueness  unless  their  distribution  and  strengths 
were  specified.)  For  homogeneous,  Isotropic  media,  (1.4)  has  been  used  to 
derive  restrictions  on  the  singularities  in  E  and  H  that  can  occur  at  an  edge. 

The  latter  were  then  used  to  prove  the  uniqueness  of  the  solution  [3].  In  this 
paper  (1. 4)  is  used  directly  in  the  formulation  of  the  problem  and  in  the  existence 
and  uniqueness  theorems. 

The  SUver-Mullo:  radiation  condition  has  been  used  to  obtain  uniqueness 
in  the  steady-state  di^actlon  problem  for  bounded  obstacles  in  a  homogeneous, 
isotropic  medium  [5,  6,8].  It  is  a  condition  which  guarantees  that  the  solution 
behaves  like  an  "outgoing  wave"  at  large  distances  from  the  obstacle.  The  form 
of  the  condition  depends  strongly  on  the  form  of  Maxwell's  equations  for 
homogeneous,  isotropic  media  and  no  such  conditions  are  known  for  inhomogeneous 
or  anisotropio  media.  However,  the  authc^  has  shown  that  for  homogeneous, 


where 
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isotropic  dissipative  media  (i,  e. ,  media  with  a  positive  electrical  conductivity) 
the  Silver-Muller  radiation  condition  implies  that  the  field  components  tend  to 
zero  exponentially  at  Infinity  [8,  p,  120].  In  particular,  it  follows  that 

for  R  sufficiently  large.  This  is  meaningful  for  Inhomogeneous,  anisotropic 
media.  Moreover,  it  is  plausible  that  in  dissipative  media  there  is  an  energy 
balance  between  the  energy  Introduced  by  the  source  fields  J  and  K  and  the 
energy  dissipated  in  the  medium,  so  that  the  time-average  energy  is  finite. 

Conditions  (1. 4)  and  (1.  5)  are  used  as  "edge  condition"  and  "condition 
at  infinity"  below.  They  can  be  combined  conveniently  into  the  single  condition 
that  the  total  ^time-average)  energy  in  the  medium  is  finite: 

(1.6)  . 

Notice  that 

(1. 7)  E,  \  =  (eJ  6  i  Ej^XEj^  -  1 E^)  =  (eJ  E^  +  E^^eJ)  +  ltE,X  '  eJ  E^) 

so  that  (1. 6)  can  also  be  written 
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§2.  Formulation  of  the  Diffraction  Problem.  The  tensors  e^j^,  and 
are  assumed  to  have  the  following  properties. 

(2.1)  €jj^(x),  |ji^j^{x)  and  are  bounded,  Lebesgue-rileasurable 

functions  of  x  c  n,  and 

(2.2)  ejj^(x),  p.^j^(x)  and  <rjjj(x)  are  uniformly  positive  definite  in  0:  1.  e. , 

there  exist  positive  constants  c  ,  u.  and  cr  such  that 

m  m  m 


fca:  all  x  <  12  and  all  real  (||1^  "  ^1  ^2  ^3^  * 

Conditions  (2.1)  and  (2.2)  imply  that  there  exist  finite  constants  c  , 

M 

UI--  and  a-.,  such 

M  M 

(2.3)  {j  5,  <  .  J^{x)  <  .r^(.  j^lx)  S,  {,) 


for  all  X  c  12  and  all  real  .  The  parts  of  conditions  (2.1)  and  (2.2)  applying 

to  €  ,  and  p, .  imply  that  the  conditions  (2,2)  and  (2.3)  on  cr.,  are 
JJC  jK  jJC 

equivalent  to  the  conditions 

However,  the  first  form  proves  to  be  more  convenient. 

The  tenscx's  and  (t^j^  are  not  assumed  to  be  symmetric.  However, 

(2.1),  (2.2)  and  (2.3)  imply  that  their  antisymmetric  parts  are  bounded  relative 
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to  their  symmetric  parts;  i.  e. ,  there  exist  finite  constants  a  such  that 


(2.4) 


l(<  ®  ^k  *  V 

^  ®  ^k  +  ’’j  V 


for  all  X  c  n  and  all  real  and  t)^  .  All  the  results  given  below  require, 
for  their  proofs,  that  (2.4)  should  hold  with  sufficiently  small  values  of  a  , 
namely 


(2.5) 


a  < 


2  V 


jn 

M 


The  formulation  of  the  diffraction  problem  given  below  makes  use  of  several 
classes  of  vector  fields  on  (2  .  To  define  them  let 


I‘,(*2)  ={A:  A{x)  is  Lebesgue-measurable  on  fl,  /  |A(x)P  dx<  «} 

^  R 

denote  the  Lebesgue  class  of  square-integrable,  complex-valued  vector  fields  on 

1  z 

a  .  Here  A(x)  =  A-(x)  +  lA  (x)  has  complex-valued  components 
Aj(x)  =  aJ{x)  +  lA^(x)  and 


|A{x)|^  =  |A\x)  +  iA^(x)i^={A\x)  +  iA^x)X-(A\x)-lA^(x))=  |a\x)|^+  |a\x)|^ 


8- 
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Notice  that,  because  of  (2.1),  (2.2)  and  (2. 3),  a  (Lebesgue-measurable) 
electromagnetic  field  has  finite  enorgy  If  and  only  if  E  c  L2(n)  and  H  c  • 

The  class  of  vector  fields  A  c  for  which  VX  A  c  12(12)  is  needed  in  the 

formulation  of  the  diffraction  problem  given  below.  Its  definition  is  motivated 
by  the  identity 

(2.6)  /a  •  VX  ®dx -/«  •  VX  Adx  =/ N  X  A  .  ®dS  , 

n  n  an 

which  is  valid  if  A  and  9  are  continuously  differentiable  and  912  is  sufficiently 
smooth. 

Defin4t^on .  Let  A  c  L2(12)  .  Then  VX  A  exists  and  equals  B  c  12(12)  4=$ 

/a*VX  *dx  =  /b  •  *dx  for  all  «c  C*(12)  . 

12  n  ° 

Op 

Here  0^(12)  denotes  the  class  of  vector  fields  on  12  which  have  continuous 

derivatives  of  all  orders  and  vanish  outside  a  compact  subset  of  12  .  VXA  is 

00 

unique,  if  it  exists,  because  (12)  is  dense  in  notations 

=  {As  A  and  VX  A  are  in  12(12)} 


and 

LJ(Vx;n)=!L2(Vx;12)r^{A:  Jk  •  VXBdx=  /b  •  VxAdx  for  all  Be  L2(Vx;12)} 

12  12  ^ 
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are  used  below.  Notice  that  "A  «  {Vx;n)"  generalizes  the  boundary  condition 

£» 

"NXA  =  0  on  an"  .  Indeed,  if  A  and  VxA  are  continuous  in  the  closure 
of  n  and  an  is  smooth  then  kt  L2(Vx;n)  Implies 

/  N  X  A  •  «  dS  -  0 

an 

for  all  4  which  are  continuously  differentiable  in  the  closure  of  n  ,  and  it 
follows  that  N  X  A  -■  0  on  an  . 

A  formulation  of  the  diffraction  problem  which  is  applicable  to  arbitrary 
domains  n,  and  media  satisfying  (2.1)  and  (2.2},  is  contained  in  the 

Definition .  Fields  E  and  H  define  a  strict  solution  of  the  steady-state 
diffraction  problem  for  the  domain  n  and  source  fields  J  e  L^Cn)  and 
K«  L2(n}^=s^£  c  L^CVxjn),  He  L2{Vx;n)  and  Maxwell’s  equations  (1.2)  hold 
almost  everywhere  in  n  . 

Notice  that  the  generalized  boundary  condition,  together  with  the  combined 
"edge  condition"  and  "condition  at  infinity"  (1.  6),  are  contained  in  the  definitions 
of  L^{'7x.',a)  and  L2(Vx;S2)  . 


10- 
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§  3.  The  Energy  Inequality.  The  fields  defined  by 


(3.1) 


K.(VXE)^+1uh.,^H^  , 


are  in  for  every  E*  and  He  L^(Vx;n),  by  (2.1).  Thus  every 

pair  E  «  L^CVxjQ),  H  «  L^CVx;^)  defines  a  strict  solution  of  the  diffraction 
problem,  and  the  correspondence  E,  H  -*  J,  K  defines  a  linear  operator  on 
L^Cn)  X  with  domain  L^CVx;  n)  X  L^CVxjQ)  .  The  main  theorem  in  this 

paper  is  an  "a  priori"  estimate  which  implies  that  this  operator  is  bounded.  It 
will  be  called 

Theorem  1  (The  Energy  Inequality).  Let  w(#0)  be  a  real  number  and  let 

*jk’  ^jk  “^jk  (2.1),  (2.2),  (2.4)  and  (2,5)  .  Then  there 

exists  a  constant  C,  depending  on  u  and  the  bounds  for  *  ,  (i,.  and  cr,, 

jk  jk  jk 

only,  such  that 


(3.2) 


Re/(*.^E.E^  +  jx,j^H.Hj^dx<  C  /(Ij1^+  lK[^)dx 


for  all  E«  L2(Vx;12),  Ht  L2(Vx;i2),  with  J  and  K  defined  by  (3.1)  .  Indeed, 
(3.2)  holds  with 


(3.3) 


C  =  4  Max(lA^  1/hiJ 


2  .  2 
<r  +  4w 
m 


m  "^m'  2  2  ,2  2. 

cj  (or  -  8watr.,  -  4a  cr. .) 

'  m  M  M' 


It  is  shown  below  that  this  number  is  positive  when  (2,  5)  holds.  Theorem  1 


and  (2,2)  imply 
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Corollarv  1.  Under  the  same  hypotheses, 


(3.4) 


/(1e[^'+  lHl^)dx<C  /(1j[^+  lK[^)dx 
S2  n 


where  =  C  Max(l/€  1/fjL^)  and  C  is  given  by  (3,3)  . 

The  proof  of  Theorem  1  is  based  on  two  lemmas  concerning  the  bilinear  form 


(3.5)  I(E,H)  =  /(E  J  -  H  K  )  dx  , 

which  may  be  stated  as  follows. 

Lemma  1.  Under  the  hypotheses  of  Theorem  1,  there  exists  a  positive  constant 

m,  depending  on  «  and  the  bounds  for  €  .  ,  n  ,  and  o-  ,  only,  such  that 

jJC  JK  JK 


(3.6) 


mRe  /(«3j,E^Ej,  +Hi^^H^H^dx<  Il(E,H)[ 


for  all  £«  L^CVxjfi),  L^CVx;^)  .  Indeed,  (3.6)  holds  with 
(3.7) 


2  „|  t  .22 

,  ,0-  -  8  w  acr  -  4a  cr 

2  2  m  M  M 

m  =  «  - - 


2  .  2 
(T  +  4« 
m 


Lemma  2.  Under  the  same  hypotheses 


|l(E,H)l  <  2Max(l/Vt^,  l/V^^)(Re /c^^E^  Ek  +  p,j^  +  WW- 


(3.8) 
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Lemmas  1  and  2  Imply  Theorem  1  with  C  =  4Max(l/€^,  which, 

with  (3,7),  gives  (3.3)  . 

Proof  of  Lemma  1.  Substituting  (3. 1)  in  (3.  b}  gives 


I(E,H)  =  /{Ej(VXH)j-H^(VXE)^+l<o«^j^E^E^-,r^j^E^lj^-i(.)jij^H^H,^}dx  . 
(2 


The  sum  of  the  first  two  Integrals  vanishes  because  E  e  L2(Vx;n),  H  i  L2(Vx;n)  . 
Thus 


UE,H)  =  \  Hj  }d*  . 

•  • 


If 


P+  —  j.  p^  p^  p”  -  p^  p^  p^  p^ 

Bji_  “  Lj  L,.  T  tij  ti.j  tji.  -  tij  —  tij  , 


‘'jk  “j  “k  "  “j  “k»  ‘'jk  “j 


=  “iic  =' “K - ' 


then  (see  (1.7)) 


(3.9)  EjE^=Ej^tiEj;.  Hj  H,  =  h;^  t  I  Hj-,  , 


and  therefore 


I(E,H)  =  +  «Pjj^H^j^dx-i/(v^j^Ejj^- wCjj^Ejj^  +  H^j^) dx 


f 


Dropping  the  underlined  terms  gives  the  Inequality 
Il(E,H)l  >(/v^j^E|^dx)2  +  2«/«^j^E"  dx/v  eJ.  dx  +  2a)  / 


+  20)  /t^jkHj^dx/o-j^Ej-j^  dx 


(3.10) 


-2“>,kE"k'^*/‘  ik^k^^  +  2a)/v^^E;^dx/^^^Ht^dx  +  a)\  /*  jj^e|^  -  K^kH^kd*)^  . 


jk“jk-^'  jk^'jk'^  "  jk'jk'^  J  '*jk"jk'- 


Now,  by  (2.4) 


’'iic‘^j  4  'fcii  =  '('jk-'w)*/  4'  -= 'jki^  ^  *  4  4>  ' 


l.e. , 


'*jk4)=i -  'fk  4k  -  »  -  “  ‘jk  4ic 'jk  4k  i  ^  -jk  4k  • 

When  multiplied  by  u  and  Integrated  over  U  this  gives 


-  a  M  /  4k  ^  <-  “  /  -  jk  4k  ^  »  M  /  -jk  4k  ^  ■ 
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Moreover,  (2. 4}  Implies  exactly  similar  results  for  and  .  Combining 
these  with  (3. 10}  gives 

W  W  ••  ee 

Eliminating  o-^j^  by  means  of  (2. 2}  and  (2.  3)  then  gives 

(3.U)  Ike.  Hll^  >,^(/.  -  4 14a.„/. 

Q  w  M 

Notice  that,  by  (3.9), 

h  »)  -  /('Jlc  ^  I'jk  • 

Hence  (3.6)  is  equivalent  to 

Thus  if 

j  then  (3.6)  follows  from  (3.11)  and  the  inequality 

I 

(3.12)  -  4l«I  a<r  j^a(af+p)  +  w^(e- P)^  >  m^(a  +  p)^ 

I 


t 
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The  proof  of  Lemma  1  will  be  completed  by  showing  that  (3,12)  holds  for  all 
real  a  and  p  when  m  is  given  by  (3.7)  .  Indeed,  (3.12)  is  equivalent  to 


(o-^-2  lul  6+  -  m^)a'  -2(  [ul  6  +  +  m^)aP  +  (a)^-m^)P^  >  0,  6  =  2a(r 


2  2  2  2  2 


M 


This  is  true  for  all  a  and  P  if  and  only  if 


(3.13) 


2  2  2 
(T  -2|wl6  +  w  -m  >0 
m  — 


and  the  discriminant 


{<r^  -  2  Iwj  6  +  co^  -  m^)  (<i>^  -  m^)  -  (  [wl  6  +  +  m^}^  >  0 


2  .  2.2 


The  last  inequality  is  equivalent  to 


(3.14)  (o- 4lwl  6  -  6^)w^  -  (o-^+ 4w^)  m^  >  0  . 


This  has  positive  solutions  m  provided 


f(6)  a  <r^  -  4Im1  6  -  5^  =  <r^  +  4«^  -  (6+ 2l«l)^  >  0  ; 
'  '  m  m 


i.  e. ,  provided  6  lies  between  0  and  the  larger  root  of  f(6)  =  0  : 

2 

0  <  6  =  2a<r  <  -2l«I+  V<r^+4w^  =  “7==*==t= 

V  O’  +  4«  +  2  Ici) 

m  ‘ 


This  inequality  is  equivalent  to  (2.5)  .  Assuming  it  holds,  (3.14)  has  solutions 

2  2  %  "  ^ 

0<m  <w  - r — . . .  6  =  2a(r., 

(T  +  4«  ^ 
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2 

Conversely,  (3.14)  Implies  (3,12)  provided  (3.13)  holds.  Now,  if  m  has 
its  largest  allowable  value  (3. 7)  , 


2  2  2 
<r  -2[c>)(6  +  c»)  -m  = 
m 


i — > 


cr  +  4w 
m 


because 


m 


6  < 


2  ^2 
cr  +  2co 
m 


+4<i)^  +  2[(o| 

m 


This  completes  the  proof  of  Lemma  1  . 

Proof  of  Lemma  2.  Notice  that,  by  (2, 2)  , 

m  m 

Thus 

IKE,H)|  =  I/(eJ.  -H.KJdxl  <  /  1e  J  ldx+  /1h  K  idx 
jj  J  1  1  J  fi  ^  Si  ‘  ^ 


<  /lElllldxt  /iHllKldx  <  (/lEl^dx)*  (/|j|^dx)^+ (/lHl^dx)^{/|Kl^dJ!)* 

n  n  Si  Si  SI  Si 

<  X(;.,,E;-^t,,^H+^dx)\/ljl^  iKl'dxji,  ^(/.j,E;t^^,H;,dx)i(/lj|^|K|W 

m  *^ni 

^■^M=’<‘''^-‘/'^X/*ik<k*>‘ilcH]kdx)M/llPt  iKl'dx)^  , 


This  is  equivalent  to  (3.8),  and  proves  Lemma  2. 
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§4.  The  Fundamental  Theorems.  Let  E,  H  and  E',  H'  be  strict  solutions 
of  the  diffraction  problem  corresponding  to  source  fields  J,  K  and  J',  K‘  , 
respectively.  Then  the  differences  E  -  E',  H  -  H'  define  a  strict  solution  with 
source  fields  J  -  J‘,  K-K',  because  Maxwell's  equations  (1.2)  are  linear. 
Hence,  Corollary  1  implies 

Corollary  2.  Under  the  hypotheses  of  Theorem  1, 

(4.1)  1Ie-e'II^+ IIh- Cj(llj-J’[|^+ (|k-k'1I^)  , 

where  (I  •  •  *  If  denotes  the  norm. 

Corollary  2  asserts  the  continuous  dependence  of  strict  solutions  on  their 
"data",  the  source  fields  J,  K  in  .  An  immediate  consequence  is 

Corollary  3  (The  Uniqueness  Theorem).  The  diffraction  problem  has  at  most 
one  strict  solution  corresponding  to  data  J  and  K  in  . 

Indeed,  if  E,  H  and  E‘,  H'  are  strict  solutions  with  the  same  data  J,  K 
then  (4. 1)  with  J'  =  J,  K'  =  K  implies  E*  =  E,  H'  =  H  .  Corollary  2  also  plays 
a  key  role  in  the  proof  of 

3 

Theorem  2  fThe  Existence  Theorem).  Let  O  be  an  arbitrary  domain  in  R  , 
let  « (#  0)  be  a  real  number,  and  let  e  p  ,  and  <r,,  satisfy  (2.1),  (2.2), 

)K  jjC  jJC 

(2.4)  and  (2.  5).  Then  the  corresponding  steady-state  diffraction  problem  has  a 
(unique)  strict  solution  for  every  pair  of  source  fields  J  and  K  in  • 

The  proof  makes  use  of  the  following  (apparently)  weaker  notion  of  solution. 
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Definition .  Fields  E,  H  define  a  solution  with  finite  energy  of  the 


steady-state  diffraction  problem  with  source  fields  J  and  K  in  9=*^  E  and 

H  are  in  (1.  e. ,  the  electromagnetic  field  has  finite  energy)  and 


(4.2)  { 


/{H.(VX  «).  -  E.  4^  -  E,  -  J.  ®.}dx  =  0  , 


f  {E  (VX  If).  +  iojLi,  .  H.  'J',  -  K.  «f.}  dx  =  0 
^  J  J  k  ) 


for  every  L^CVx;!})  and  'S' e  . 

The  usefulness  of  this  notion  stems  from 

Lemma  3.  Fields  E,  H  define  a  solution  with  finite  energy  E,  H  define 
a  strict  solution. 

Proof.  (The  implication  "  -srr-p  ")  Note  that  C^{n)C.L^{Vy.;i2)C.L^('7x;a)  . 
Hence,  identities  (4. 2)  imply 


and 


/e  •  Vx  iS'dx  =  /(-iujA  H  +  K  I’S  dx 

n  JK  K  J  J 


-09 


for  all  *  and  *  in  C.(n),  where  the  fields  iu*,,  E,  +0-,,  E,  +  J  and 

O'  ”  jk  k  jk  k 

+  Kj  are  in  L^Cn)  .  This  Implies  (i)  VXH  and  VXE.  exist  in  L^Cn) 
and  (il)  Maxwell's  equations  (1.2)  hold  almost  everywhere  in  Q  .  In  particular 
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E  and  H  are  in  L^CVx;  12)  .  Finally,  the  second  of  the  identities  (4.2), 
with  Kj  =  (VxE)j  +•  iwjjijj^  Hj^,  gives 

/ {E,(VX  'if).  -  (VXE).  )dx  =  0  for  all  St  c  L2(Vx;  12)  ; 

12 

i.  e. ,  E«  L^(Vx:12)  .  Thus  E,  H  defines  a  strict  solution. 

(The  implication  "  <=  ")  Multiply  the  first  of  Maxwell's  equations  (1.  2)  by 
4  €  L^(Vx;12)  and  integrate  over  12  .  Then  the  term 

/(VXH)  4  dx  =  /h  (VX  4)  dx 
12  12  •'  ■' 

because  4  t  L2(Vx:  12)  and  H  e  L2(Vx:  12)  .  This  gives  the  first  of  Identities 
(4. 2)  and  the  second  follows  by  a  similar  argument. 

Maxwell's  equations  (1.2)  determine  a  linear  operator  A  on  the  Hilbert 
space  L2(n)  X  L2(0),  with  domain 

D(A)  =  L®(Vx;12)  X  L2(Vx;12)  , 

defined  by  A  (A,  B)  a  (  J,  K)  where 

c 
I 
{ 

! 

i 


Corollary  1  implies  that  A  is  one-to-one  on  D(A),  and  that  A  ^  is 

2.  Z  2t 

bounded:  I|A,  b1|  <  Cj  1|a(A,  B)|| ,  where  i|A,  b1I  =  UaH  +  UbU  is  the  norm  in 
L2(n)XL2(fi)  . 

Theorem  2  is  equivalent  to  the  assertion  that  R(A),  the  range  of  A,  is 

equal  to  X  •  This  will  be  established  by  showing  that  (i)  R(A)  is 

dense  in  X  L2(^)  (li)  R(A)  is  closed  in  X  •  The  first 

i' 

assertion  is  called 

Theorem  3«  Under  the  hypotheses  of  Theorem  2, 

^  X  , 

where  the  bar  denotes  closure  in  L^CO)  X  . 

Proof.  R(A)  is  a  closed  linear  subspace  and 

L^{a)  X  =  RUT  ®  N  , 

by  a  standard  theorem  on  Hilbert  space  [2,  p.  25]  .  Hence,  it  is  sufficient  to 
show  that  if  E,  H  €  R(A}  is  orthogonal  to  R(A)  in  L^CSi)  X  L^Cn)  then 
E  =  H  =  0  .  Now  E,  H  _L  R(A)  means 

/{E  f  +  H.  K.)dx  =  0  for  all  A  c  L°(Vx;n),  B  c  L,(Vx;n)  . 

J  J  J  1  ^  ^ 

Combine  this  identity  with  Maxwell's  equations  (4.  3),  and  take  first 
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and  second 

A  =  0,  B  =  *  with  »  t  L^CVx; «)  . 

This  gives  the  pair  of  identities 

+  H^(VX  *)^}dx  =  0  , 

/{E.(VX  «-).  -  >I'^}dx  =  0  , 

for  all  ««L2(Vx;Q),  «  L^CVx;^)  .  These  identities  state  that  E,  H  Isa 

solution  with  finite  energy  of  a  modified  diffraction  problem  (the  adjoint  problem), 
with  -<o  for  u,  source  fields  J  =  K=0  . 

But,  the  conditions  for  the  validity  of  Theorem  1  clearly  imply  the  same  conditions 
for  the  modified  problem.  Thus,  E,  H  is  a  strict  solution  of  the  modified  problem 
(by  Lemma  3)  with  source  fields  J  =  K  =  0  and  therefore  E  =  H  =  0  by  Theorem  1. 
This  proves  Theorem  3. 

Proof  of  Thecrem  2.  Let  J  t  L^CO),  K  e  L^CO)  .  Then,  by  Theorem  3,  there 
exist  sequences  of  fields  e"  €  L^{Vx;Q),  h"  t  L^CVx;!!)  whose  source  fields 
j”,  K*'  converge  to  J,  K  in  L^CfJ)  X  L^Cn)  .  Applying  Corollary  2  with  E  =  e”  , 

H  =  h",  E*  =  e"*,  H*  =  H*"  gives 

He"  -  E'"i^+  Dh"  -  <  Cj((1j”  -  ff  +  |k"  -  k“1I^)  . 

It  follows  that  {e’^}  and  {h’^}  define  Cauchy  sequences  in  L2(i2}  •  Hence, 


limit  fields 
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E  =  lim  e",  H  =  lim  (in  L^CO)) 

n-*«o  n  -*«o 

n 

exist  because  L2(i2)  is  complete.  Now,  E  and  H  define  strict  solutions 
with  source  fields  .  Hence,  by  Lemma  3,  E*^  and  define  solutions 

with  finite  energy;  i. e. ,  identities  (4.2)  hold  with  E*'’  for  E,  h”  for  H  , 
j”  for  J  and  k’’  for  K  .  Making  n-*’*  in  these  identities  gives  the  same 
identities  for  the  limit  fields,  since  all  the  fields  converge  in  L-(i2)  .  Thus 

£0 

E,  H  is  a  solution  with  finite  energy,  and  therefore  a  strict  solution  (by  Lemma  3) 
having  the  prescribed  source  fields  J,  K  which  proves  Theorem  2. 
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